Draft version February 1, 2008 

Preprint typeset using I£T^X style cmulatcapj v. 7/15/03 



O 
O 

(N 

5-1 . 

< 

ON 



Or 
i < 

o ■ 

to: 



> 



^1" 

O 

o 



IDENTIFICATION OF ABSORPTION FEATURES IN AN EXTRASOLAR PLANET ATMOSPHERE 

T. Barman 

Lowell Observatory, 1400 W. Mars Hill Rd., Flagstaff, AZ 86001, barmanaiowell.edu 

Draft version February 1, 2008 

ABSTRACT 

Water absorption is identified in the atmosphere of HD209458b by comparing models for the planet's 
transmitted s pectrum to rec e nt, m ulti- wavelength, eclipse-depth measurements (from 0.3 to 1 fim) 
published by iKnutson et al.1 (|2007l ). A cloud-free model which includes solar abundances, rainout 
of condensates, and photoionization of sodium and potassium is in good agreement with the entire 
set of eclipse-depth measurements from the ultraviolet to near-infrared. Constraints are placed on 
condensate removal by gravitational settling, t he bulk metallicity, a nd th e redistribution of absorbed 
stellar flux. Comparisons are also made to the ICharbonneau et al.l (|2002f) sodium measurements. 
Subject headings: planetary atmospheres - extrasolar planets 



1. INTRODUCTION 

The discovery of transiting extrasolar planets has 
opened the door to direct detections and characteriza- 
tion of their atmospheres. Observations using the STIS 
instrument on HST provided the first glimpse of what 
the photospheric composi tion is like for a nearby EGP. 
ICharbonneau et "all (|2002h measured the relative change 
in eclipse depth for HD209458b across a sodium doublet 
(5893A) resulting in the first detection of atomic absorp- 
tion in an EGP atmosphere. Following the Na detec- 
tion, IVidal-Madjar et"all (|2003h discovered an extended 
hydrogen-rich atmosphere surrounding HD209458b using 
a similar technique as Charbonneau et al., but in the UV. 
At Lyman-a wavelengths, HD209458b is ~ 3 times larger 
than in the optical. These detections were made using a 
technique called transit spectroscopy as the planet passes 
in front of its star. Transit spectroscopy uses the fact that 
the wavelength-dependent opacities in the planet's atmo- 
sphere obscure stellar light at different planet radii lead- 
ing to a wavelength-dependent depth of the light-curve 
during primary eclipse. Consequently, searching for rela- 
tive changes in eclipse depth as a function of wavelength 
directly probes the absorption properties of the planet's 
atmosphere with the potential to reveal the presence (or 
absence) of specific chemical species. 

In this paper, recent measurements of HD209458b's 
radius are combined in a multi-wavelength comparison to 
model atmosphere predictions. Atmospheric molecular 
and atomic absorption are identified and constraints are 
placed on the basic atmospheric properties. 

2. THE LIMB MODEL 

Transit spectroscopy probes the limb of a planet which 
is the transition region between the day and night sides. 
One would, therefore, expect that, in the presence of 
a horizontal temperature gradient between the heated 
and non-heated hemispheres, the temperatures across 
the limb w ould be cooler t han the average daysi de tem- 
peratures ([Barman et alj l2005t llro et all 120051 ). Re- 
cent Spitzer observations of both transiting and non- 
transiting hot-Jupiters showing large flux variations with 
phase have provided strong evidenc e supporting such a 
day-to-night temperature gradient (ICharbonneau et al.l 
l2005UDeming et al.ll2005bl 120061: lHarrington et aH 2006). 



Describing the limb ultimately requires a multi- 
dimensional model atmosphere solution; however, as is 
common practice, a simpler one-dimensional model is 
used here to represent the average properties of the limb, 
in both a longitudinal and latitudinal sense. To explore 
a variety of limb temperature structures, the incident 
stellar flux has been scaled by a parameter a. A model 
with a = 0.25 represents an average description of the 
entire planet in t he presence of very e fficient horizontal 
energy transport (|Barman et alj I2005D . Increasing a to 
0.5 increases the heating of the model atmosphere mak- 
ing it more appropriate for just the dayside which, in this 
work, will represent an upper limit to the plausible mean 
temperatures at the limb. 

The basic chemistry across the limb is modeled assum- 
ing chemical equilibrium, determined by minimizing the 
free energy while including grain formation. To under- 
stand the effects of gravitational settling of grains on the 
chemistry and the transmission spectrum, the removal 
of grains from the atmosphere is included via two simple 
approximations that represent opposite e xtremes. The 
first is the "cond" approximat ion used by iBarman et al.l 
(|2001[) and lAllard et al.l (|200lD which simply ignores the 
grain opacity without altering the chemistry or abun- 
dances. The second is the "rainout" approximation 
which iteratively reduces, at each layer, the elemental 
abundances (by the appropriate stoichiometric ratios) 
involved in grain formation and recomputes the chem- 
ical equilibrium with each new set of stratified elemen- 
tal abundan ces. This approach is similar to other rain- 
out m odels (|Feglev fc Loddersl Il996t iBurrows fc Shard 
1999), except that the depletion of elements is contin- 
ued until grains are no longer present. 

The transmission of stellar fluxes through the limb of 
HD209458b is determined by solving the spherically sym- 
metric radiative transfer equation while fully accounting 
for scattering an d absorption of both intrinsic and ex- 
trinsi c radiation |Hauschildt] Il992t IBarman et alj 1200 lL 
2002). Spherical, instead of the more traditional plane- 
parallel geometry, naturally accounts for the curvature of 
the atmospheric layers and changes in chemistry along 
the slant paths through the upper atmosphere. The 
planet radius at a given wavelength (Ra) is obtained by 
determining the radial depth at which the transmitted 
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Fig. 1. — Monochromatic transit radii over the STIS spectral range for the baseline rainout model with (red) and without (blue) 
photoionization of Na and K. The solid and dotted red lines are the same, except H2O line opacity is excluded in the latter. Horizontal 
bars correspond to mean radii across bins with A-ranges indicated by the width of each bar. S TIS measurements by Knutson et al. (2007) 
are shown in green with 1 a error bars. Vertical dashed lines mark the narrow A-range used by Charbonncau ct al. (2002). 



flux is equal to e 1 times the incident starlight along 
that same path. 

3. RESULTS 

A cloud-free atmosphere with rainout, a = 0.25, and 
solar abundances is adopted as the baseline model. This 
model, along with o thers, is compa r ed to the relative 
Ra measurements of iKnutson et al.l (|2007f ) which have 
a reported precision high enough to constrain many of 
the basic atmospheric properties. Since a comparison is 
being made to relative Ra values, the model results were 
uniformly scaled to match the observations in the 4580 
to 5120A wavelength bin; this scaling was always less 
than 0.005R,/ U)9 . Overall, the baseline model (red solid 
line in Fig. [T]) reproduces the observed rise in Ra toward 
shorter wavelengths, the increase across the Na doublet, 
and the increase at the far red wavelengths. The baseline 
model comparison to the data has a x 2 that is 3 times 
smaller than a constant Ra- 

3.1. Water Absorption 

Water is predicted to be one of the most abundant 
species in an EGP atmosphere and, given its broad ab- 
sorption features in the infrared, plays a crucial role in 
regulating the temperature-pressure (T-P) profile. The 
first major H2O absorption band appears between 0.8 
and 1 f jm, a region covered by the last two wavelength 
bins of IKnutson et al.l (|2007f ) . As illustrated in Fig. QJ 
there is excellent agreement between the baseline model 
and the observations in this part of the spectrum espe- 
cially across the longest wavelength bin that sits on top 
of the H2O band. Qua litative l y simi lar w ater features are 
seen i n the models of iBrownl ()200lD and iHubbard et aT] 
(|2001h : however these models fall many a below the ob- 
servations. The baseline model also predicts mean Ra- 
peaks equal to 1.330, 1.343, and 1.341 Rj up for the next 
three water bands (at A ~ 1.15, 1.4, and 1.9 yum). 

A model that excludes H2O line opacity is also shown 
in Fig. Q]and is greater than 10 a below both the observa- 
tions and the baseline model prediction. Removing H2O 



opacity also produces a significant drop in Ra near 0.9 
/jm that further increases the discrepancy between the 
model and overall red/near-IR observations. No other 
opacity source could be responsible for the observed rise 
in Ra across this part of the spectrum. 

3.2. Photoionization 

After the repo rted sodium detection in th e atmosphere 
of HD209458b (jCharbonneau et all l2002f) . there were 
several attempts to explain why the strength of this fea- 
ture was much lower than expected based on earlier mod- 
els (e .g., by Seager & Sasselov, 2000). iBarman et al.1 
(2002) explored departures from local thermodynamic 
equilibrium (LTE) which are capable of producing an in- 
version in the cores of the N a doublet line profiles . How- 
ever, the earlier models of IBarman et al.l (|2001l I2002D 
were constructed under the cond approximation and, 
consequently, contained a larger number of free met- 
als along with TiO and VO molecular absorption com- 
pared to a rainout model. These additional sources of 
optical/UV opacity lead to a hotter upper atmosphere 
and a very shallow photoionization depth for Na (only 
a 5% reduction of Na across the regio n probed by tran- 
sit spectroscopy) . iFortnev et al.l (|2003l ) also explored Na 
photoionization and found that their atmosphere model 
(which included rainout) could be brought into reason- 
able agreement with the Na observations. The models 
presented here account for the angular dependence of 
ionziation on the limb's day s ide, b ut not on the night 
side. However. IFortnev et al.1 (|2003l ) have shown that Na 
ionization is still present at ~ 5° past the terminator for 
pressures relevant to the transit spectrum modeled here. 

Figure [T] compares cloud- free solar abundance rainout 
models with and without photoionization of Na and K. 
While these models include photoionization, the LTE ap- 
proximation on the atomic level populations and line 
source function is maintained. Ionization (radially) at 
the limb reaches 50% at P ~ 0.1 and 1.2 mbar for Na 
and K, respectively, and stops at P ~ 2 and 7 mbar, re- 
spectively, resulting in a reduction of the RA-peaks across 
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Fig. 2. — Relative flux differences (with time from center of 
eclipse) between a wavelength band centered on the Na doublet and 
the mean of two wavelength bands on either side of the Na doublet. 
The total wavelength range is indicated by vertical lines in Fig. [T] 
Each model is cloud-free with solar abundances and a = 0.25. See 
Figure legend and text for the distinguishin g characteristics of each 
model. Symbols are the observations of Charbonncau ct al. (2002) 
with 1 <t error bars. 



the Na and K lines. For K, this reduction extends out to 
the line-wings resulting in a significantly smaller mean 
Ra, bringing the model into ~ 3a agreement with the 
Knutson et al. measurement. The impact of Na pho- 
toionization is mostly confined to the core of the doublet 
leading to only a small reduction of the mean Ra , but suf- 
ficient to bring the model into ~ la agreement with the 
observations at these wavelengths. Though not included 
here, ionization past the terminator onto the night side 
should further improve the agreement across the Na and 
K doublets. The two narrow features on the red wing 
of the K doublet are due to Rb, which should also be 
affected by photoionization due to its very low first ion- 
ization potential. Photoionization of Rb resulted in a 
near complete removal of these features, but reduced the 
mean Ra for this bin by less than 0.002 Kj up - 

The broad wavelength bins allowed Knutson et al. to 
obtain very precise relative Ra measurements; however, 
as illustrated by the Na doublet, such broad bins lim- 
its the constraints that can be placed on atomic absorp- 
tio n features. In contrast, th e narrow range analyzed 
by ICharbonneau et aTl (|2002D resulted in much larger 
Ra error-bars, but still precise enough to easily distin- 
guish between the various models shown in Fig. [TJ The 
flux differences across the Na doublet as a function of 
time (measured from the transit center) were computed 
for each model using the sam e narrow wavelength bins 
as ICharbonneau et aTl (|2002l ). indicated in Fig. [1] by 
vertical dashed lines. Fig. [2] compares model transit 
curves to the 2002 Na measurements and shows a large 
discrepancy between the observations and solar abun- 
dance models with pure equilibrium chemistry (cond or 
rainout). Including photoionization brings the baseline 
rainout model into rough agreement with the observa- 
tions. Note the cond model shown in Fig. [3] does not 
include photoionizat ion (similar to the LTE model from 
I Barman et al.l (|2002f )) and illustrates how far off the pre- 
dicted Ra can be across the Na doublet under simplified 
assumptions. 

3.3. Rainout, Metallicity, and Temperature 



The rainout model used here reduces the individual 
metal abundances with depth iteratively until clouds do 
not form, thus mimicking efficient gravitational settling. 
However, the removal of metals is not 100% efficient, 
leaving behind a variety of atoms (in addition to Na and 
K) to contribute to the line opacity. The impact of rain- 
out on the planet's atmosphere is made apparent by com- 
paring rainout and cond models (Fig. [3]). In the cond 
model grain opacities are simply ignored in the transfer 
equation while the chemistry remains that of a pure equi- 
librium model without any actual removal of refractory 
elements. Thus, the cond model represents the minimum 
impact of grain formation on the stratified abundances 
and leaves behind a considerable amount of free metals 
along with molecules like TiO and VO. This leads to 
much stronger atomic lines, including Na and K along 
with TiO and VO bands. 

Apart from a very minor contribution from molecules, 
the Ra features at A < 0.8/im are all due to atomic line 
opacity 1 . These narrow features are due to blended lines 
from metals like Ca, Al, Fe, Ni, Mn, and Cr. Note that, 
while Rayleigh scattering does contribute to the opacity 
in the blue/UV, a removal of atomic line opacity would 
drop Ra by ~ 0.02R,/ UJ , for A < 0.45/im. The top panel of 
Fig. [3] compares the solar abundance rainout model with 
photoionization (shown in Fig. [IJ to a rainout model 
with 10 x solar abundances. Since the metal abundances 
were uniformly scaled, the grain formation and corre- 
sponding removal of refractory elements were also uni- 
formly enhanced. The additional opacity altered the T-P 
profile as well as the total atmospheric extension. These 
factors contribute to a similar Ra pattern in the opti- 
cal but enhanced molecular absorption features in the 
red/near-IR. The increase in metals increases the H2O 
absorption feature along with the wings of the Na and 
K doublets (which form deeper in the atmosphere where 
photoionization is less important). Also, FeH and CrH 
absorption is larger in the metal-rich atmosphere leading 
to an increase in Ra on both sides of the K doublet. The 
metal-enhanced Ra spectrum is in good agreement with 
the observations near 0.9/im, but is noticeably too high 
across both the K doublet and the H2O band. 

Transit spectroscopy can also help constrain the tem- 
peratures across the limb. The lower panel of Fig. [3] 
compares an a = 0.5 model (i.e. a model with a hotter 
dayside-like T-P profile) to the cooler a = 0.25 base- 
line model. In the hotter model, grain formation is less 
pronounced resulting in a greater concentration of most 
free metals and, thus, stronger UV/blue absorption lines 
than in the a = 0.25 model. Higher temperatures at 
depth also result in equilibrium concentrations of Na 
and K that are several times lower than found at the 
cooler temperatures of the a = 0.25 model. The net re- 
sults are two distinctively different Ra spectra with the 
a = 0.25 model being more consistent with the Knutson 
et al. measurements. 



4. SUMMARY 



1 IBallester etUI ( f2007T) claim the R A rise in the blue/UV is due 
to the Blamer edge of hot hydrogen in an extended atmosphere. 
While the present model does not support this explanation, if ver- 
ified, it could indicate that rainout is more efficient than predicted 
here. 
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Fig. 3. — Monochromatic transit radii over the STIS spectral range. Top panel: solar abundance models, cond (black) and rainout(red), 
are compared to a model with 10 X solar metal abundances including rainout (blue). Lower panel: solar abundances baseline rainout model 
with a = 0.25 (red) compared to a similar model with a = 0.5 (blue). Horizontal bars have the same meaning as in Fig. [T] 



Photoionziation plays an important role for both Na 
and K, and potentially many other species. No evidence 
is found to support a large metal enhancement (e.g., 10 x 
solar), though smaller Jupiter- like enhancements are not 
ruled out. Furthermore, the agreement between model 
and observations demonstrated here alleviates the need 
for substantial cloud coverage along the limb between 
0.05 and 0.001 bars, whic h would other wise truncate and 
flatten the K\ spectrum (|BrownH2001l ). In the baseline 
model, the predicted location of clouds (e.g., MgSiOs) 
lies just below the minimum Ra(^ l-315Rj up ) across the 
STIS wavelength range, consequently deep clouds (and 
also very high clouds) remain a possibility. 

The models presented here also predict R> variations 
(~ 0.02 Rj up ) across the 2-0 R-branch of CO in the K- 
band. This is inconsistent (at 2.5a) with a non-detection 
of CO b ased on Keck- NIRSPEC transit spectra taken 
in 2002 (|Deming et alJ l2005a[ ) . It is likely that using a 
single T-P profile to represent the horizontal and pole- 
to-pole variations across the limb averages out too much 
of the upper atmospheric structure; this simplification 
might explain the K-h&nd discrepancy. High clouds may 
also be involved. In addition, the Keck and HST obser- 
vations were taken between 1 and 2 years apart and time- 
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dependent at mospheric variation s along the limb cannot 
be ruled out (iMenou et al.ll2003f) . 

While iKnutson et alJ <|200tT) did not attribute their 
measured Ra variations to any absorption features (this 
was not the focus of their paper), the models presented 
above clearly show that these measurements are consis- 
tent with strong water absorption near 1 /im. A detection 
of water in the limb of HD209458b is nominally at odds 
with a recent Spitzer I RS spectrum that shows no H2O 
features for this planet dRichardson et al.ll2007D . These 
data were taken during secondary eclipse and directly 
probe the planet's dayside with negligible contribution 
from the limb. It is possible that the day side atmosphere 
is nearly isothermal (jFortnev et al.ll2006l ) which would re- 
sult in a spectrum with no detectable water absorption 
features, despite a copious water supply. The transmis- 
sion spectrum, however, would contain absorption fea- 
tures independent of an isothermal dayside or limb. 
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